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' Phase transitions of a flux lattice in layered superconductors with magnetic field perpendicular 

Q\ , to the layers and in presence of disorder are studied. We find that the Josephson coupling between 

Q\ ' layers leads to a strongly pinned Josephson glass (JG) phase at low temperatures and fields. The JG 

t-H | phase undergoes either a decoupling transition with increasing field or a depinning transition with 

increasing temperature. The resulting phases undergo further depinning and decoupling transitions, 
respectively, resulting in a phase diagram with a multicritical point where four phases meet. The 
, phase diagram accounts for unusual data on B^SV^CaCi^Os such as the "second peak" transition 

and the recently observed depinning transitions. 
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The phase diagram of layered superconductors in a magnetic field B perpendicular to the layers is of considerable 
interest in view of recent experiments on high temperature superconductors [lj. A first order transition in Y Ba2Gu^O-j 
q ■ (YBCO) and in BiiSriC aCuiO% (BSCCO) has been interpreted as a melting transition of the flux lattice. This 
first order transition terminates at a multicritical point, which for BSCCO is at Bq sa 300 — 10 3 G and To w 
P ( l 40 - 50A, while for YBCO f| it is at B « 2 - 10T and T « 60 - 80A, depending on disorder and oxygen 
concentration. The multicritical point also terminates a "second peak" transition which is manifested by a sharp 
increase in magnetization; the transition line at B sa Bq and T < Tq is weakly T dependent. Neutron scattering 
and fiSK data |D,|5| show that positional correlations of the flux lattice are significantly reduced near these phase 
boundaries, except however, near the multicritical point where a reentrant behavior is observed ||. Recent data on 
Ndi.s5Ceo.i5Cu04s (NCCO) has also shown a second peak transition; here, however, B decreases with temperature 
i-^j | near the superconducting transition at T c « 23 A with no apparent multicritical point. 

In a recent remarkable experiment Fuchs et al. ^] have shown that the phase diagram of BSSCO is much more 
elaborate. They show that the spatial distribution of an external current exhibits a transition from bulk pinning 
to surface pinning of vortices with most of the current flowing at the sample edges. This depinning line crosses the 
multicritical point and its temperature is almost B independent at B < Bq . The depinning transition correlates with 
7-H ! anomalies in vibrating reed experiments (^] and in magnetization [jl0|| . Thus there are four phase transition lines 
^ ■ which emanate from the multicritical point at Bo, Tq\ The first order line, the second peak line and depinning lines 
for both B < B and B > B . 

An extensive theoretical effort has been devoted to understanding the field-temperature (B — T) phase diagram |ll]] 
in presence of disorder. In particular, it was proposed that at low T and B a Bragg glass is stable fl^Jl3|l , exhibiting 
algebraic decay of translational order and divergent Bragg peaks Jl4| ). Melting is expected to occur by thermal or 
j^. ' disorder induced dislocations, as indeed demonstrated for fields parallel to the layers j^,[u|. 

^\ , The flux lattice can undergo a transition which is unique to layered superconductors, i.e. a decoupling transition 
In this transition the Josephson coupling between layers vanishes while the lattice can be maintained by 
the electro-magnetic (e-m) coupling between layers. A disorder induced decoupling was also proposed as a crossover 
phenomena [T^ j. 

The theory of Daemen et al. jL8) employed the method of self consistent harmonic approximation (SCHA) to find 
the decoupling temperature Td{B). The SCHA leads to a conceptual difficulty since it predicts that the Josephson 
coupling vanishes for all purposes at T > T^. Koshelev [^0| has shown that above some critical temperature the 
£j , Josephson critical current vanishes, however, a finite Josephson coupling is maintained and in fact accounts for 
the experimentally observed plasma resonance. Thus the decoupling transition, as found by SCHA, needs to be 
reinterpreted. 

In the present work we consider temperatures below the melting temperature T m of the flux lattice and study (i) 
the decoupling transition in a renormalization group (RG) framework and (ii) effects of disorder by employing replica 
symmetry breaking (RSB) methods. We find a glass transition T g such that for T < T g strong pinning is expected. 
The lines Td and T g cross and lead to four distinct phases which meet at one point in the B — T phase diagram, 
remarkably close to the experimental phase diagram [ j]-|4|,|s|,|lC| . 

Consider a flux lattice with an equilibrium position of the Z-th flux line at R/ . The flux line is composed of a sequence 
of singular points, or "pancake" vortices, whose positions at the n-th layer can fluctuate to R/ + u". Consider the 
transverse part of u™ with the Fourier transform ux{(\,k), where q, k are wavevectors parallel and perpendicular to 
the layers, respectively. The elastic energy due to e-m coupling is given by 
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H e - m = lY^ida 2 ) 2 ^ 2 + clMkl\\u T {^k)\ 2 (1) 

q.A: 

where the flux line density is 1/a 2 , d is the spacing between layers, q is within the Brillouin zone (of area (27r/a) 2 ), 
fc| < n/d and k z = (2/d) sin(fcd/2). The shear and tilt moduli are given (for a ^> d) by J21J 

c° 6 = r/(16da 2 ) 

4M = [r/(8da 2 X 2 ab k 2 )} ln(l + a 2 k 2 J^) (2) 

where r = ip^d / (4tt 2 X 2 b ) sets the energy scale and A a (, is the magnetic penetration length parallel to the layers; 
r w fO 3 - 10 4 K for YBCO or BSCCO parameters §. Note the strong dispersion of c% A {k) which decreases by the 
large factor (d/a) 2 when k varies from k < 1/a to 1/a < k < ir/d. 

The Josephson phase between the layers n and n + 1 at position r in the layer involves contributions from a 
nonsingular component and from singular vortex terms. The singular phase around a pancake vortex at position 
R; + u" is a(r — R; — uf) where a(r) = arctan(y/a;) with r = (x, y). We assume that all vortices belong to the flux 
lines, i.e. there are no free pancake antipancake {pp) pairs which appear as relevant fluctuations only above T m . The 
effect of the nonsingular component is a negligible T/r term in the RG equation p2| while expansion of the interlayer 
phase difference a(r — R; u") — a(r — R; u" +1 ) yields for the Josephson phase b n (r) — X^( u ™ +1 — u")Va(r — R;). 
The Hamiltonian is then 

n = n e - m - (j/^o) Yl I d2r cos 5 «( r ) ( 3 ) 

n 

where J is the interlayer Josephson coupling and ^0 is the coherence length, serving as a short distance cutoff. Since 
Vq ~ 1/r decays slowly, even if u™ are small the contribution of many vortices which move in phase (q — > 0) 
leads to a divergent response of 6„(r). In Fourier space, the relevant &(q, k) fluctuations involve q < 1/a where 
6(q, k) = 2nid{e lkd — l)uT(q, k)/q, i.e. enhanced q — > fluctuations. 

Standard RG proceeds |^2| by integrating high q components leading to a new cutoff £ > £0 and a £ dependent 
coupling </(£). The significant softening of (fc) at k > 1/a implies that the k integration is dominated by k w ir/d 
so that the resulting Cg 6 [(7 4 /fc 2 ]|6(q, k)\ 2 term from Eq. (1) can be replaced by an upper cutoff on the q integration, 
q u = 21og 1/2 (a/d)/A a6 . To first order in JJT we obtain J(£) - {£q u ) [1/2{1 ~ t)] where t = T/T d and the decoupling 
temperature (similar to the SCHA result ]lq] ) is 

_ 4a 4 /■ dfc x ^ ra 2 log(a/d) 
U ~ cP [ J cUk)> ~ 4.A 2 b (4j 

Thus for T > Td </(£) vanishes on long scales (£ — > 00). Second order RG results in renormalization of C44 and in 
generation of Josephson coupling between next nearest neighbors ]22j ] . The second order terms enhance Td by a factor 
(1 — jJ/T)^ 1 where 7 is a nonuniversal parameter. 

The RG process shows that the decoupling transition is manifested only on long scales. Thus the thermal average 
of the local observable < cos6„(r) > remains finite at T > Td, as in the J/T expansion p0| . This, however, does not 
imply long range order in cos6„(r) - the same J/T expansion yields a power law decay for < cos 6„(r) cos 6„(0) > 
correlation, as also obtained from RG. 

A hallmark feature of two-dimensional superconductivity is the In p dependence of a pp interaction on their sepa- 
ration p, leading to a power law I-V relation ]2^ j. To probe this feature in the decoupled phase we consider a high 
temperature expansion of Eq. (R|) with an added pp pair. The direct e-m interaction is [|22f ~ hip while the additional 
effective free energy to order J is 

F p¥ (p)~(J 2 /T) fd 2 r [ d 2 /|r-rr 4t [l-cosK(r,p)-a (r',p))] (5) 

J J\r-r'\>(o 

where cto(r 7 p) — a(r — p) — a(r). Eq. ([5]) can be shown to be bounded by a ~ In 2 p term, supporting a nonlinear 
I-V relation at T > Td- In contrast, at T < Td the pp interaction increases as ~ p leading to a finite critical current. 
Thus the decoupling transition is manifested by the change in correlation function, vanishing of the Josephson critical 



current |20| and by nonlinear I-V relation. 

We proceed now to study effects of disorder. Since T < T m we assume first small fluctuations |u™| <C a. Consider 
a short range pinning potential U™ in (r) which couples to the vortex shape function p(r) as J d 2 r^2 n z U™ in {v)p{r — 
R; — u") . Expansion in u™ and averaging U™ in (r) by the replica method p3] leads to the replicated Hamiltonian, 
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H T jT=\ V [ C (fc)gX/3-s &6 Q (q,fc)6 /3 *(q,fc) 

ft/ — 

q,fe;a,/3 2 



-7^2 £ / ^cos6«(r)-p ^ / d 2 rcosK(r)-^(r)] (6) 

where a, (3 are replica indices, c(k) = (a 2 /2Trd) 2 c1 4 (k)/T and So — Ua 2 d/ (Aird 2 T) 2 with [7 an average of the pinning 
potential. In Eq. (^) the Cg 6 term has been replaced by a cutoff q u on q integrations, as above. The inter-replica 
Josephson coupling, i.e. the v term in Eq. (||), is generated from the J term in second order RG. It is essential to keep 
the v term from the start since it couples different replica indices and can lead to distinct physics by RSB jl6],^5) . 

Note that the more general form of the disorder term is Jl2|-jl4| cos[Q • (u"' a — u"'^)] where Q is a reciprocal lattice 
vector; expansion of this cosine leads to the so term in Eq. (^)7^ The cosine form is essential for deriving the Bragg 
glass properties of the flux lattice, i.e. the 1/r decay of the displacement correlation at distances r > £. The domain 
size £, over which the flux lattice is well correlated will be of significance below. 

The RSB method Q proceeds by employing a variational free energy T va r = J~o+ < H — Ho > with To the 
free energy of Ho — \ J2q k-a p G~^(q, fc)6"(q, /c)6 /3 *(q, k) and G Q ^(q, k) is determined by an extremum condition on 
F var . This yields 

k ) = I c (% 2 + Z 1<W " so{q 2 /k 2 z ) - a a ,p (7a) 
z = (J/2T£ 2 d) exp[-§ ]T G a a (q, k)} (7b) 

<r a ,0 = (v /£,od)[exp(-±B at/3 ) - S a ^ ^ ex -P(~^ B a n )} (7c) 
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where B a ^ = 2^ q fc[G a;Q (q, k) — G Qj /3(q, k)] and z is a renormalized Josephson coupling. The method of RSB |24|] 
represents a hierarchy of matrices such as <7 Qj/ 3, B Q ,/3 in terms of functions <j(u), B(u), respectively, with < u < 1. 
The amount by which the replica symmetry is broken is measured by a glass order parameter A(«) = ucr(u) — 
<t{v) dv. Using standard methods [^,^5) we find that the solution for A(u) is a step function, i.e. A(it) = for 
u < 2t while A(u) = Ao for 2t < u < 1, where 

(z + AoVA^^/^A^ 1 - 2 *) (8) 



with the cutoff A c « c{ir/d)q 2 L . Thus a solution with A ^ is possible only if t < 1/2. To solve for z in Eq. (7b) we 
need the diagonal part, 

G a , a (q, k) = \n{2etv/ze a d) + (s /8tt 2 )[I(z) + zl'(z)] (9) 

where I(z) = j d 2 qdk/{k 2 c(k)[c(k)q 2 + z]} and I'(z) ~ dl(z)/dz. Formally I(z) diverges at k = 0; this divergence 
can be traced back to our assumption that the cos[Q • (u"' Q — u™ ,/3 )] term is expanded into the sq term in Eq. (|^). 
Retaining this cosine leads to Imry-Ma type domains of correlated u™ whose size perpendicular to the layers is £ z . 
Within a domain the u™ expansion is valid so that ir/£ z serves as a lower cutoff in the k integration. More formally, 



keeping the cos[Q • (u™' a — u™ )] term replaces so in Eq. (7a) by a matrix a a p which corresponds to an RSB function 
<7i(u). This leads to an additional glass order parameter Ai(u) = uai(u) — J Q " <t\(v) dv and the divergent I(z) is 
replaced by a term in Eq. (0) of the form 



1 fUv &l(u)q 2 /k 2 z 



(i k c(k)q 2 +z J v 2 c(k)q 2 + z + A (it) + A 1 (u)q 2 /k 2 



which converges at k — » 0. The general solution for both A(u) and Ajjw) involves a rather difficult set of two coupled 
differential equations. For J = v = the Bragg glass solution is Jl^ , p^ | Ai(u) ~ u 2 for u < u c and Ai(w) = Ai(w c ) 
for u c < u < 1, where u c ~ 1/so is small. Thus for t = T/T^ ^> u c the structure of Ai(u) at small u should not be 
affected by A(it) with its step at u — 2t. The scale at which the k divergence is cutoff is at k < Ai(it c )/c(0) ks \ jl z . 

Consider then I(z) with a lower cutoff tt/£ z on the k integral. If l z < a it leads to a small correction ~ 0(d/£ z ) to 
the main 1/a < k < n/d integration range. If £ z > a then the ir/£ z < k < 1/a range in I(z) can be neglected if 

£ z < (a 4 /d 3 )/[32/n 2 (a/d)]. (11) 
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Since o>dwe expect this to be valid; in particular for BSCCO parameters (see below) it is always valid. Note also 
that a finite thickness of the sample can also serve as a cutoff replacing £ z . 

I(z) with l/a<fc<7r/d integration yields (so/8ir 2 )I(z) = 2s log(A c /z) where the dimensionless disorder parameter 
is s = AirU\\ h /[T 2 a 2 log 2 (a/d)] . The renormalized Josephson coupling of Eq. (7b) is then 



z/A c = e- 1 [J 2 /(8T 2 *eo^A c )] 1 /( 1 - 2s ). (12) 

Comparing Eqs. (|s|]l2|) shows that Ao vanishes at s = t (up to a nonuniversal ~ 1/hxv term) and formally there is a 
solution with Ao < when s < t. However, the average distribution |m| of o(q, k)\ 2 is ~ exp[— |6(q, k)\ 2 / G a ^ a (q, k)] 
is acceptable only if G a ^ a {(\,k) > 0. This is a thermodynamic stability criterion and for our solution it reduces to 
Ao > 0. Thus the regime where both z, Aq are finite is limited to s < i, t < s; we term this regime the Josephson 
Glass (JG) phase. The glass parameter vanishes (continuously) at t — s while the Josephson coupling vanishes 
(continuously) at s = |. For s > | and t < \ the solution is z — while Ao ^ satisfies Eq. (||), i.e. it is a decoupled 
glass phase. Finally, for Ao = 0a replica symmetric solution is valid at s < t < 1 — s with 

z/A c » ( J/2T^ 2 rfA c ) 1 /( 1 -^*) . (13) 

Thus s + 1 = 1 for s < | defines a decoupling transition. 
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FIG. 1. Phase diagram. Full lines are decoupling lines where the Josephson coupling vanishes. Dashed lines are depinning 
lines where the Josephson glass parameter vanishes. B is determined by the disorder strength while T = \Td(a = ao) (Eq. (Q)) 
where Oq = Bq/^q. 



The phase diagram, Fig. 1, has four phases which all meet at a point defined as Bq,Tq. Bq is determined by the 
disorder strength via s = \ while T = \Td{a = a ) (Eq. (|j)) where a 2 , = B /<j>o. Since s increases with B the 
s — i line defines a decoupling transition from a JG phase at low B to a pinned glass phase (G) at high fields. It is 
remarkable that although the G phase has vanishing Josephson coupling (z = 0) the Josephson induced disorder (the 
v term in Eq. (^|)) is dominant in determining the glass nature of this phase. In fact, RG shows that J first increases 
(scaling from £o to l/q u ), generating the v term, and only at scales beyond l/q u J decreases to zero. 

The JG phase at B < Bo undergoes another transition at t = s, i.e. at T — Tq (up to Yn.B factors) into a phase with 
finite Josephson coupling while the glass parameter Ao vanishes. This Josephson (J) phase has, however, the Bragg 
glass type disorder. The condition that t z has negligible effects in the JG or G phases implies that the pinning effect 
from the Josephson induced Ao is much stronger than that associated with the Bragg glass. Thus the JG-J transition 
is a depinning transition, from strong to weak pinning. The G phase also undergoes a depinning transition into a 
decoupled Bragg glass phase (D) at T = BqTq/B; the D phase is a Bragg glass phase maintained by the interlayer 
e-m coupling. 

The J phase undergoes a decoupling transition at B — 2BqTq/(T + Tq), using s w B/2Bq. The J-D transition is 
continuous, at least for small J/T where RG and RSB are valid. The SCHA has been formally extended to higher 
J/T and found to be of first order fig] . We recall now that melting and related dislocations have been neglected. An 
upper bound on T m is the vortex transition p2| at r/8 where pp pairs are thermaly excited. Thus To < r/8 limits 
our description near Tq to ao ^ X a b- 
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We reconsider now the condition for the Bragg glass with its domain size l z to have a negligible effect on our 
glass phase. The coefficient of the cos[Q • (u™' Q — u™ ,/3 )] term is related to so and yields the Bragg glass length E^] 
i z « 0.004A^ fc /[sdlog 2 (a/d)]. This domain size has non-negligible effects only if l z > a and also if t z does not satisfy 
Eq. (|l|), i.e. 2 ■ I0 3 ald/\l b < a < 0.5d\ ab /a with s w a„/a 2 . This range exists if a < 0.06A ab , i.e. B > 10 5 G for 
A a6 w 17001. Thus Bragg glass effects are negligible for BSCCO and YBCO, except near B in the B w 10T YBCO 
sample 0. Note also that even if ag < 0.06A a b, Bragg glass effects can be neglected if the sample thickness satisfies 
Eq-©- 

We interpret the experimentally observed second peak phenomena [pHwl as the JG-G transition, i.e. a decoupling 
transition within the glass phase. While a decoupling scenario has been suggested as a crossover phenomena [p],^9|, 
the present theory predicts a strict phase transition. The JG-G transition at B = Bq is T independent up to Tq and 
Bo decreases with impurity strength; both features are consistent with experimental data JjJ-0]. The NCCO sample H] 
has similiar parameters to BSCCO samples implying a comparable To; however, the low value of T c k 23K indicates 
that a multicritical point with Tq < T c is probably not realized. 

Recent data [^]-[l0) has shown an additional phase boundary in BSCCO, i.e. a depinning transition line which 
crosses the critical point Bo, Tq. Our result for the depinning temperature, T = To at B < Bq being B independent 
(up to ~ lni? terms), is in accord with the data. At B > Bq we expect the depinning line at T = BqTq/B, in 
qualitative agreement with a stronger B dependence. [^0[ . 

Neutron data |6) has shown a reentrant behavior in the 600 — 10 3 G range with positional correlations increasing 
with temperature. This is consistent with our decoupled phase which is weakly pinned, leading to enhanced positional 
correlations. The reentrant behavior seems to extend to B < Bo so that our J-D line may be the first order line, at 
least near Bq; at lower fields this decoupling line probably joins the melting line. 

In conclusion, we have found a phase diagram which is remarkably close to the experimental one fl|41Jsl,|i"o| , having 
a multicritical point where four phases meet. Our theory provides a fundamental interpretation of both the second 
peak transition and the recently discovered depinning transition. 
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